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Abstract 

It has been shown recently that a simple formula in terms of the valence nucleon numbers and the 
mass number can describe the essential trends of excitation energies of the first 2 + states in even- 
even nuclei. By evaluating the first order energy shift due to the zero-range residual interaction, 
we find that the factor which reflects the effective particle number participating in the interaction 
from the Fermi orbit governs the main dependence of the first 2 + excitation energy on the valence 
nucleon numbers. 
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Valence nucleoli numbers N p and N n have been extensively adopted to parameterize var- 
ious nuclear properties phenomenologically. The valence proton (neutron) number N p (N n ) 
is defined as the number of proton (neutron) particles below the mid-shell or the number of 
proton (neutron) holes past the mid-shell within the given major shell. Hamamoto was the 
first to show that the square root of the ratio of the measured and the single-particle B(E2) 
values, [B{E2) exp /B{E2) sp] ^ j is roughly proportional to the product N p N n Casten 
extended the idea and suggested the N p N n scheme where he showed that if we parameterize 
the collective variables or operators in even-even nuclei in terms of the product N p N n then 
we get a substantial reduction in the number of parameters without serious loss of accu- 
racy . It was also demonstrated that the N p N n scheme could be applied not only to 
even-even nuclei but also to odd- A and odd-odd nuclei {4]. 

Recently, we proposed a simple empirical formula in terms of the valence nucleon numbers 
N p and N n , and the mass number A for the excitation energies E x (2f) of the first 2 + states 
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FIG. 1: Excitation energies of the first 2 + states in even-even nuclei. The solid lines represent 
isotopic chains. Part (a) shows the data quoted from Ref. 6, and part (b) shows the results 
obtained by Eq. (TTJ . 
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in even-even nuclei jsj. The formula is given by 

E x = aA'i + (3 [exp(-Ai\y + exp(-AJV„)] (1) 

where the parameters a, (3, 7, and A are fitted from the data. From Fig.[TJ which is repro- 
duced from Ref. 5, we find that the essential trends of the excitation energies E x (2f) are 
well reproduced throughout the periodic table. The data displayed in part (a) of Fig.[T]is 
quoted from the recent compilation by Raman et al. |6|] and the results shown in part (b) 
of Fig.ODwere obtained by applying Eq. ([!)) where a = 34.9 MeV, (3 = 1.00 MeV, 7 = 1.19, 
and A = 0.36 [5|. 

In this brief report, we want to examine what makes such an odd dependence of the 
excitation energies E x (2f) as given by Eq. ([TJ) on the valence nucleon numbers N p and N n . 
Since we are dealing with the systematics observed throughout the whole periodic table 
instead of inspecting the detailed structure of a single nucleus, we want to keep the physics 
involved as simple as possible. For that purpose, we first find the difference between the 
energy of the lowest two quasi-particle state and the excitation energy of the first 2 + state. 
Then we compare it with the first order energy shift due to the residual interaction. Here, 
we obtain the former energy difference from the measured quantities and the latter energy 
shift we evaluate by using simple perturbation theory. 

Within the standard BCS theory [3], the single quasi-particle energy E k of the k-th orbit 
is given by 



E k = v/^-z^ + A 2 (2) 

where e k , fi, and A are the single particle energy, the Fermi energy, and the gap energy, 
respectively. We can obtain the gap energy A empirically from the binding energies of the 
adjacent nuclei. Let B(Z, N) be the negative binding energy of a nucleus whose atomic 
number and neutron number are Z and N, respectively. Then the proton gap energy A p 
and the neutron gap energy A n can be approximated by the following three point expressions 



A p « l - \B{Z + 1,N) — 2B(Z, N) + B(Z-1,N)\, (3) 
A n w l - \B(Z, N + l)- 2B(Z, N) + B(Z, N - 1) | . (4) 

Let us write the excitation energy E x (2f) of the first 2 + state by simple perturbation 
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FIG. 2: Upper part shows twice the gap energy 2A, which is the lesser of 2A p (Eq. ([3])) and 2A n 
(Eq. (HI)). Binding energies are taken from Ref. 10. Lower part shows the measured excitation 
energies of the first 2 + states (the same as those shown in Fig.QJa)). 

theory as 

E x {2t)=E^+E^ + --.. (5) 

Then the unperturbed energy E^ is equal to the energy of the lowest two quasi-particle 
state, Ek + Ey where k and k' indicate the orbits which belong to the lowest two quasi- 
particle state. Also the first order energy shift E^ is given by 

E {1) =< 3k3k'\V\j k j y >j. =2 + . (6) 

which is the expectation value of the residual interaction V on the lowest two quasi-particle 
state which is coupled to the total angular momentum J w = 2 + . 

From Eq. ([2]), we have that E k ~ A for the single particle orbits k near the Fermi level. 
Therefore, we take twice the gap energy 2 A in place of the lowest two quasi-particle energy 
E k + -Etc'191]. At the upper part of Fig. [2, twice the gap energy, 2 A, is shown for even-even 
nuclei. Here 2A is determined by the lesser of 2A p and 2A n which is estimated empirically 
either by Eq. ([3]) or by Eq. [lOj . We can observe that all of the 2 A are located higher 
in energy than the measured E x (2f) data, which is shown at the lower part of Fig.[2j We, 
therefore, want to compare the difference between 2 A at the upper part and E x (2f) at the 
lower part of Fig.[2]with the first order energy shift E^ which will be given shortly. 

In order to estimate E^\ we take only the central part of the residual interaction in the 
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zero-range approximation, 



V(l,2) = -v 6 3 (r 1 -r 2 ), 



(7) 



where vq, which is the strength of the zero-range interaction, will be treated as the only 
model parameter in our calculations. Analytic expressions for the matrix elements of the 
zero-range interactions are readily available from the literature llj]. The expectation value 
appearing in Eq. (jf)]) is given by 



< jkjk'\V\jkjk' >J-=2+ 

<R 4 >(2j k + l)(2j k , + l] 



-Vq 



1 1, T 9 

< Jk-Jk' - -| JO > 



407T (1 + S jkjk> ) 
x {[1 + (-) ifc+4 '](«, + «,) + [2 - (-l) jk+jk '](u Al + v A ,f] 



+ (-) jk+j *' < jk-jk'-\Ji > 2 (u jk v jk , + (-; 



(8) 



where < R 4 >, v ? , and v3 are the integral of the radial wave functions 

' Jk '» 



< i? 4 >= 



drr z 



R h( r ) R h> 



(9) 



the occupation probability of the fc-th orbit, and 



respectively. 



We evaluate the matrix element in Eq. ([SJ) by using the following procedure. First, the 
single-particle wave functions and energies are generated by diagonalizing the Woods-Saxon 
(WS) potential 

Vws(r) = ^r^r (10) 



1 + exp 

where for the WS parameters we use Vq = 50MeV, Rq = 1.27fm, and a = 0.67fm. Then 
the radial integral < R A > of Eq. (jHJ) is evaluated using the WS wave functions. Also the 
occupation probabilities u| fc in Eq. (JSJ) are determined by solving the BCS equations with 
the simple monopole pairing interaction 12| 

1 



< (jij2)J\v pah \(j 3 j 4 ) J >= --GS Jfi ^(2j 1 + l)(2j 3 + l) 



(11) 



where the strength G of the pairing interaction is adjusted to reproduce the measured gap 
energy by Eq. or Eq. (pE} for each nucleus in our calculations. 

In Fig.[3], we compare our calculated results (solid squares) for —E^ with the measured 
differences (solid circles) 2 A — E x (2f) for the following four cases: (a) Z = 52 isotopic 
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FIG. 3: Calculated results (solid squares) for the energy shift are compared with the mea- 

sured differences (solid circles) 2A — E x (2f) for the following four cases: (a) Z = 52 isotopic chain, 
(b) Z = 70 isotopic chain, (c) N = 70 isotonic chain, and (d) N = 102 isotonic chain. 

chain, (b) Z = 70 isotopic chain, (c) N = 70 isotonic chain, and (d) N = 102 isotonic chain. 
In our calculations, the strength w of the zero-range interaction in Eq. ([7]) is adjusted so 
that the calculated —E^ coincides exactly with the measured difference 2 A — E x (2±) for 
the heaviest lead isotope ^fPb in Fig. [31(d). The resulting adopted value of v in all of our 
subsequent calculations for E^ 1 ' is vq = 1, 700 MeVfm 3 . By observing the results displayed 
in Fig. [31 we find that the first order energy shift —E^ calculated by such a simple model 
like perturbation theory describes the measured difference 2 A — E x (2f) unbelievably well. 
Thus, this results enable us to examine E w in place of E x (2f) to find the origin of its 
dependence on the valence nucleon numbers N p and N n . 

If the total angular momentum jf of the Fermi orbit is larger than 1/2, which holds 
practically in most cases of our calculations, then both of the two quasi-particles which form 
the lowest two quasi-particle state belong to the Fermi orbit. Therefore jk and jy in Eq. (jSJ) 
can be replaced by jf. This, in turn, reduces the expression for the energy shift m 
Eq. flE]) to a very simple form which is written as 



E {1) = <jfjf\V\jfjf -./•: 2- 



20tt < >< 3f 2 Jf 



12, >' 



where we use the known value of the Clebsch-Gordan coefficient 
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FIG. 4: Calculated results (solid squares) for the energy shift -E^ are compared with the effective 
particle number (2jf + l) 2 u| u| for the same four cases as those in Fig.0 Note that the scale for 
the energy shift is shown at the far left axis while that for the effective particle number is shown 
at the far right axis. 

which holds identically for any even J and half integer values of j. 

The expression for the first order energy shift can be divided into two major factors 
as shown by the large square brackets in Eq. (II 2 1) . One is the radial integral < i? 4 > times the 



Clebsch-Gordan coefficient squared and the other is (2j /+l) 2 u 2 t>| which can be interpreted 
to reflect the number of effective particles that participate in the interaction from the Fermi 
orbit. Out of the two factors, we expect that, the former is kept more or less the same 
over different isotopes and isotones. Therefore, the variation in values of follows that 
of the effective particle number (2jf + l) 2 u 2 f v 2 f . In order to check this point, we compare 
the negative of the first order energy shift —E^ (solid squares) with the effective particle 
number (2jf + l) 2 u? f v 2 (solid circles) in Fig.H] for the same four cases as those in Fig.0 
Note that the scale for the energy shift is shown at the far left axis while that for the effective 
particle number is shown at the far right axis in Fig.Hl We can confirm undoubtedly from 
Fig. HI that the behavior of E^ follows that of the effective particle number quite closely. 

Finally, we compare the measured excitation energies (solid circles) E x (2f) of the first 2 + 
states in Fig. [5] for the same four CclSGS clS those in Fig. [3] with the one calculated by perturba- 
tion theory (solid squares), namely 2A + E^\ and with the one obtained from the empirical 
formula (solid triangles) given by Eq. ([T]). By observing these plots, we can ascertain, first of 
all, that the empirical formula, Eq. ([TJ, expressed in terms of the valence nucleon numbers 
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FIG. 5: Measured excitation energies (solid circles) E x (2f) of the first 2 + states for the same four 
cases as those in Fig. [3] are compared with the one calculated by the perturbation theory (solid 
squares), namely 2A + E^\ and with the one obtained from the empirical formula (solid triangles) 
given by Eq. (TjQ). Measured excitation energies E x (2f) are quoted from Ref. 6. 

describes the measured excitation energies E x (2f) quite well, even quantitatively to some 
extent. Furthermore, we can also confirm that simple perturbation theory can be employed 
in examining the possible dependence of E x (2f) on the valence nucleon numbers at least 
qualitatively. 

In summary, we have examined the recently proposed empirical formula, Eq. (pQ), for the 
possible origin of its dependence on valence nucleon numbers. Recently, it has been shown 
that Eq. ([T]), which depends on the mass number A and the valence nucleon numbers N p 
and N n in a very simple fashion, can describe the essential trends of excitation energies 
E x (2f) of the first 2 + states in even-even nuclei throughout the whole periodic table [5]. In 
order to find out what makes such a dependence of ^(2^) on N p and N n , we calculated 
the first order energy shift resulting from the zero-range residual interaction by using 
simple perturbation theory. Then the shif was compared with the difference between the 
energy of the lowest two quasi-particle state and the measured excitation energy E x (2i). 
The lowest two quasi-particle energy was approximated by twice the gap energy determined 
from the binding energies of the adjacent nuclei through the three point expression given 
by Eq. (j3J) or (jlj). In our calculations, the strength vq of the zero-range residual interaction 
given by Eq. (J7J) was the only free parameter. It was kept fixed in all of our calculations 
at Vq = 1, 700MeVfm 3 for which the calculated —E^ coincided exactly with the measured 
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difference 2A — E x (2f) for the heaviest lead isotope 1 g|Pb. We found that the variation in 
values of followed that of the factor (2jf + l) 2 w| w| which can be interpreted as the 
effective particle number participating in the interaction from the Fermi orbit. Therefore, 
we concluded that the effective particle number governed the main dependence of the first 
2 + excitation energy on the valence nucleon numbers as given by Eq. (pQ). 
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